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ABSTRACT 


The  design  and  performance  analysis  of  a  three-channel  Robust 
Symmetrical  Number  System  (RSNS)  virtual  spacing  direction  finding  (DF) 
system  was  carried  out  in  this  thesis.  The  design  was  based  on  hardware  used 
in  the  previous  work  in  which  a  digital  three-channel  RSNS-based  interferometer 
was  built.  In  designing  a  RSNS  virtual  spacing  DF  system,  it  was  found  that  the 
ratio  between  the  modulus  is  an  important  parameter.  The  ratio  affects  the 
selection  of  the  moduli  set,  the  array  size  and  the  system’s  field  of  view.  Based 
on  the  hardware  used,  the  system’s  signal-to-noise  ratio  (SNR)  was  determined 
for  a  simulated  emitter.  Random  Gaussian  noise  was  injected  into  the  system  to 
determine  the  accuracy  of  mapping  the  actual  angle-of-arrival  (AOA)  to  the 
estimated  AOA.  Monte  Carlo  simulations  were  carried  out  to  determine  the 
probability  of  correctly  estimating  the  AOA  of  the  signal.  Simulation  results  from 
the  effect  of  spacing  errors,  phase  errors,  frequency  error,  moduli  set  change  and 
SNR  change  was  presented  and  analyzed.  The  RSNS-based  DF  system  was 
found  to  be  relatively  insensitive  to  spacing  error.  A  proposed  two  set  of  moduli 
were  used  to  improve  the  accuracy  and  probability  of  correct  estimation  in  both 
low  SNR  and  high  SNR  conditions. 
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I.  INTRODUCTION 


A.  DIRECTION  FINDING 

Radio  direction  finding  is  a  class  of  direction  finding  (DF)  by  which  the 
direction  to  a  microwave  signal  is  determined  by  means  of  a  radio  direction  finder 
(RDF).  A  RDF  is  a  receiving  system  that  analyzes  the  information  extracted  from 
the  passing  electromagnetic  radio  wave  to  obtain  direction  of  arrival  (DOA) 
information  [1].  In  military  context,  some  of  the  more  important  uses  of  direction 
finding  relating  to  electronic  warfare  applications  include  early  warning  threat 
detection,  targeting  and  electronic  intelligence  [2]. 

An  RDF  could  be  separated  into  three  broad  categories  depending  on  the 
method  by  which  it  determined  the  angle  of  arrival  (AOA)  of  the  signal.  The  three 
categories  are  [2]: 

1.  Amplitude  Comparison 

2.  Phase  Delay 

3.  Time  Delay 

Amplitude  Comparison  uses  either  a  direct  or  comparative  method  to  obtain  AOA 
information  from  the  amplitude  response  of  the  antenna  elements. 

Phase  Delay  and  Time  Delay  systems  extract  the  AOA  information  by  a 
similar  principle  except  that  the  former  used  phase  information  while  the  latter 
used  the  time  of  arrival  information  of  the  wave.  Both  systems  required  at  least 
two  separate  antennas  with  pre-determined  spacing.  Figure  1  showed  how  a 
phase  delay  or  time  delay  method  would  work.  The  plane  wave  would  arrive  at 
antenna  element  1  earlier  than  antenna  element  2.  The  antenna  spacing,  d, 
would  cause  a  time  delay  or  phase  delay  between  the  antenna  elements. 
Assuming  that  the  frequency  of  the  wave  and  the  antenna  spacing  are  known, 
the  AOA  can  be  determined. 
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Incident  Plane  Wave 


Figure  1 .  Basic  phase-delay  /  time-delay  interferometry  (After  [1]). 


To  increase  resolution,  a  larger  antenna  spacing  is  usually  used. 
However,  a  larger  antenna  spacing  also  introduces  ambiguities.  If  the  spacing  is 
greater  than  a  half-wavelength  distance,  more  than  one  unique  AOA  exists. 

B.  PREVIOUS  WORK 

Through  the  years,  research  into  using  a  Robust  Symmetrical  Number 
System  (RSNS)  for  direction  finding  has  been  investigated  in  order  to  achieve 
high  resolution  with  a  small  antenna  baseline.  It  has  been  shown  that  a  three- 
channel  DF  array  based  on  RSNS  theory  and  using  a  digital  architecture  can 
achieve  a  good  resolution  of  1 .25  degrees  RMS  with  a  baseline  of  66.2  cm  [3]. 

C.  PRINCIPAL  CONTRIBUTIONS 

This  thesis  focused  on  the  design  and  analysis  of  a  three-channel  virtual 
spacing  DF  system  with  two  levels  of  processing.  Design  considerations  for  the 
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RSNS  virtual  spacing  DF  system  was  analyzed.  It  highlighted  the  importance  of 
selecting  the  appropriate  modulus  ratio  for  sizing  of  the  antenna  array  as  well  as 
field  of  view  requirements. 

By  analyzing  the  hardware  used  in  the  previous  implementation  of  the 
RSNS-based  DF  system,  the  receiver’s  noise  figure  was  determined.  From  the 
noise  figure,  a  link  equation  was  used  to  calculate  the  theoretical  signal-to-noise 
ratio  (SNR)  of  a  simulated  emitter.  Simulations  were  carried  out  using  this  link 
equation  in  determining  the  performance  of  a  three-channel  RSNS  virtual 
spacing  DF  system.  The  effects  of  spacing  errors,  phase  errors,  frequency 
errors,  moduli  set  change  and  SNR  change  were  determined  from  the  results  of 
the  simulation.  Monte  Carlo  simulations  were  carried  out  to  determine  the 
probability  of  correctly  estimating  the  AOA  of  the  signal.  Simulation  results  from 
the  effect  of  spacing  errors,  phase  errors,  frequency  error,  moduli  set  change  and 
SNR  change  was  presented  and  analyzed.  The  RSNS-based  DF  system  was 
found  to  be  relatively  insensitive  to  spacing  error.  A  proposed  two  set  of  moduli 
were  used  to  improve  the  accuracy  and  probability  of  correct  estimation  in  both 
low  SNR  and  high  SNR  conditions. 

D.  THESIS  OUTLINE 

Chapter  II  reviews  the  basic  operating  principles  of  RDF  systems  and 
provides  the  fundamental  basis  for  applying  the  RSNS  for  direction  finding. 

Chapter  III  introduces  the  RSNS  theory  and  how  it  can  be  applied  to 
phase  interferometry.  Previous  experimental  results  are  also  presented. 

Chapter  IV  explains  how  the  RSNS  interferometry  works  with  virtual 
resolution. 

Chapter  V  presents  the  considerations  for  designing  a  three-channel 
RSNS  virtual  spacing  direction  finding  system. 

Chapter  VI  explains  the  models  used  in  the  simulation. 

Chapter  VII  presents  the  results  and  analysis  of  the  simulations. 

Chapter  VIII  presents  the  conclusion. 
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II.  PHASE  DELAYED-BASED  INTERFEROMETRY 


This  chapter  reviews  the  basic  operating  principle  of  phase  delayed-based 
interferometry  and  how  AOA  ambiguities  can  arise. 

A.  THEORY 

It  is  assumed  that  the  DF  array  operates  in  the  far  field;  hence  the  incident 
electromagnetic  wave  is  planar.  A  one-channel  linear  interferometer  is  shown  in 
Figure  2.  The  two  antenna  elements  are  space  a  distance  d  apart  and  the 
incident  plane  wave  arrives  at  an  angle  0  as  shown.  The  angle  of  arrival  can 
vary  over  a  range  of  ±90° . 


Figure  2.  Two-element  interferometry  (From  [4]). 
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The  incident  wave  arrives  at  antenna  element  1  first  and  after  an 
additional  distance  of  dsin(#),  it  arrives  at  antenna  element  2.  A  sinusoidal 
(time-harmonic)  signal  received  by  the  antenna  elements  has  the  following  form: 

V-,  =  V  cos  (a>ct  +  y/j  +<?>j)  (2.1) 

where  /  denotes  the  element  index  number;  V  is  the  signal  amplitude  response, 
a>c  is  the  signal  frequency  in  radians/sec,  y/j  is  the  phase  term  due  to  the  cables 


and  (f>j  is  the  phase  term  due  to  path  delay  relative  to  the  reference  antenna 

(antenna  1).  The  received  signals  are  fed  into  an  analog  phase  detector  which 
consists  of  a  mixer  and  a  low-pass  filter  (LPF).  After  the  phase  detection  mixing 
process,  the  output  voltage  has  the  following  form: 

i/2 

Vmix W)  =  Vl(MV2(</>2’t)  =  -^- [C0S(2 coct  +  X Y  +  X ^  +  cos( A ^  +  A^)]  (2 -2) 

The  term  A y/  results  from  the  phase  differential  due  to  cable  length  differences. 
The  term  A(f>  results  from  the  phase  differential  due  to  the  additional  d sin(<9) 
distance  traveled  by  the  signal  wave-front  to  reach  antenna  element  2  and  it  can 
be  written  as 


A </>  =  -  fa  =  kd  sin(6>) 


where  k  = 


T 


and  A  is  the  wavelength. 


(2.3) 


After  the  LPF  removes  the  high  frequency  component  of  the  signal,  the 
output  signal  which  remains  is 

V2 

voiyfW  =  -^-cos(A^  +  A^)  (2-4) 

We  can  measure  Ay/  and  eliminate  it.  After  normalizing  the  voltage’s  amplitude, 
the  final  output  signal  is 

Vout(0)  =  cos(A^)  =  cos(/rdsin(60)  (2.5) 
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B.  AMBIGUITY 


As  the  phase  differential  term  repeats  itself  with  a  period  of  ,  it  can 
result  in  an  ambiguity  if  distance  between  the  antenna  elements  is  greater  than 

half  a  wavelength.  To  illustrate  this  point,  we  first  consider  a  case  where  d  =  ^ , 

which  gives  Atf>  =  ;rsin(6>) .  There  will  be  a  one-to-one  mapping  of  the  AOA  to  the 
phase  difference  as  shown  in  Figure  3  and  no  ambiguity  exists. 


Angle  of  Arrival  (degrees) 

Figure  3.  Phase  differential  vs  AOA  for  d  =  j/^  (From  [5]). 

Next,  we  consider  a  case  where  d  =  A,  which  gives  A#)  =  2;rsin(0) .  For 
every  phase  difference  measured,  there  will  be  two  possible  AOAs  as  shown  in 
Figure  4. 
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Angle  of  Arrival  (degrees) 


Figure  4.  Phase  differential  vs  AOA  for  d=X  (From  [6]). 

The  accuracy  of  the  AOA  measurement  can  be  stated  as  the  change  in 
phase  difference  for  a  given  change  in  AOA.  Hence  from  Equation  (2.4), 

^^-  =  kdcos(6)  (2.6) 

dO 

For  any  given  frequency  and  AOA,  Equation  (2.6)  shows  that  the  accuracy  of  the 
measurement  is  directly  proportional  to  the  distance  d.  Hence,  a  longer  baseline 
would  give  a  more  accurate  measurement. 

To  overcome  the  ambiguity  coupled  with  the  need  for  higher  angular 
resolution,  dual  baselines  are  often  used.  The  longer  baseline  is  used  to  provide 
the  high  angular  resolution  required,  while  the  shorter  baseline  is  used  to  resolve 
any  ambiguities  [2]. 

C.  FOLDING  WAVEFORMS 

The  output  voltage  from  any  two  antenna  elements  can  be  represented  by 
Equation  (2.5).  For  the  cases  of  d  =  j/^  and  d  =  A ,  the  normalized  voltages  are 
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plotted  in  Figure  5.  It  shows  that  the  both  waveforms  are  symmetrical  about  the 
zero  angle  of  arrival  (broadside)  and  the  visible  region  is  from  ±90  degrees.  The 
number  of  folding  periods,  n,  within  the  visible  region  is  given  by 


d  2d 


(2.7) 


Figure  5.  Output  voltage  vs  AOA  ford  =  ^  and  d  =  A  (From  [3]). 

In  a  binary  interferometer  system  where  each  successive  baseline  is 
changed  by  a  factor  of  two,  it  is  possible  to  resolve  an  ambiguity  of  the  longer 
baseline  and  yet  achieve  higher  resolution  [2].  A  three-channel  interferometer  is 
shown  in  Figure  6.  The  successive  baseline  is  changed  by  the  factor  of  2. 
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1 


2 


3 


4 


Y  Y  Y  Y 

- - %= A - ] 

* - d}=2A - - 


Figure  6.  Three-channel  interferometer  with  binary  relationship 

between  spacing  (From  [2]). 


The  size  of  the  physical  antenna  element  array  is  often  the  limitation  in 
determining  the  accuracy  of  the  DF  system.  The  next  chapter  introduces  another 
method  of  increasing  resolution  without  the  requirement  of  a  larger  antenna 
element  baseline. 
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III.  ROBUST  SYMMETRICAL  NUMBER  SYSTEM-BASED 

DIRECTION  FINDING 


As  it  is  advantageous  to  have  a  higher  resolution  without  the  need  for  a 
large  antenna  baseline,  a  signal  preprocessing  method  utilizing  the  robust 
symmetrical  number  system  (RSNS)  was  introduced  to  DF  in  [3-7],  This  chapter 
reviews  the  RSNS  theory  and  how  it  is  implemented  in  phase  interferomentry. 

A.  THE  ROBUST  SYMMETRICAL  NUMBER  SYSTEM 

The  main  idea  behind  RSNS-based  direction  finding  is  to  represent  a 
symmetrical  waveform  out  of  an  antenna  by  a  RSNS  sequence.  The  RSNS  is 
based  on  a  sequence  defined  as 

xm  =  [0,1,2,....,m-1,m,m  +  1,....,2,1]  (3.1) 

where  xm  is  a  row  vector  and  m  is  the  modulus  and  is  a  positive  integer  (m>0). 

In  an  A/-channel  RSNS  [8],  where  N>2,  the  basic  sequence  for  the  /th 
channel  with  modulus  m  is 


The  sequence  is  repeated  in  both  directions  and  it  forms  a  sequence  with  a 
period  of 

Prsns  =  2  A/my  (3.3) 

Considering  a  single  channel,  the  discrete  states  of  the  RSNS  can  be  expressed 
as 


n-Sj 

N 


9  =  i 


2Nmj  +  N-n  +  Sj  - 1 
N 


Sj  <  n  <  Nmi  +  s;-  + 1 

Nmi  +Sj  +2<  n  <  2/Vm,-  +  s,-  - 1 


(3.4) 


where 
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g  is  the  nth  level  of  channel  /, 
rrij  is  the  /th  channel  modulus, 

S/  is  the  corresponding  sequence  shift  s,-  =  0,1,2, _ , A/  —  1  (mod  N)  and 

N  >  2  is  the  number  of  sequences  in  the  system. 

The  values  (s1,s2,....,sA/)  must  form  a  complete  residue  system  modulo  N 

and  because  of  the  relative  property  of  the  shifts,  one  of  the  shift  values  will  be 
set  equal  to  0. 

An  /V-channel  RSNS  is  formed  by  selecting  N  moduli  and  N  shift  values 
Sj ,  (1  <  /  <  N) .  Since  the  fundamental  period  for  channel  i  is  2A/n?,- ,  the  RSNS 

vectors  must  be  a  mutilple  of  2A/m/ .  Hence,  the  fundamental  period  for  the 
RSNS  is 

pfRSNS  =  {2mlN,2m2N! 2mNN}  (3.5) 

where  {a-|,a2,....,aw}  is  the  least  common  multiple  of  a-\,a2,....,aN .  From  Equation 
(3.5), 

pfRSNS  =  2N{m^,m2,...,mN)  (3.6) 


To  illustrate  the  RSNS,  consider  a  modulus  set  of  [3  4]  where 
m-)  =3(s-|  =0)  and  m2  =4(s2  =1).  Table  1  shows  the  sequences  illustrating  the 
inherent  integer  Gray-code  property. 
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Table  1.  RSNS  sequence  for  m- |  =3(s-|  =0)and  m2  =4(s2  =1) 

(From  [6]). 
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The  system  dynamic  range,  M ,  is  the  maximum  number  of  distinct  vectors 
without  any  ambiguity.  From  Table  1,  the  dynamic  range  for  [3  4]  is  M= 15.  It 
was  shown  that  the  shifts  employed  will  affect  the  dynamic  range  sequence 
starting  and  ending  points  but  it  does  not  have  any  effect  on  the  value  of  M  [9]. 

The  dynamic  range  of  several  moduli  in  combination  has  been  formulated 
in  closed  form  [10].  For  two  channels,  there  are  three  cases  of  two  relatively 

prime  moduli  and  m2  that  yield  M . 


Case  1:  For  m- |  >  3  and  m2  =  m- 1  +1 , 

M  =  3(/77-|  +  m2 )  -  6 

(3.7) 

Case  II:  For  /r?1  >  5  and  m2  =  m-\  +2, 

M  =  3(/77-|  +m2)-l 

(3.8) 

Case  III:  For  /r?1  >  5  and  m2  =  m-\  +C,  and  C  >  3, 

M  =  4  m.|  +  2  m2  -  2 

(3.9) 

For  the  example  in  Table  1,  which  satisfies  Equation  (3.7), 

M  =  3(/771  +  m2 )  -  6  =  3(3  +  4)  -  6  =  1 5 . 

For  three  channels,  there  are  two  cases  of  closed  form  solution  that  have 
been  formulated  [11,  12]. 

•  Case  I:  For  =  2k  - 1,  m2  =  2k  + 1  for  >  3, 

M  =  +  7  (3.10) 

•  Case  II:  For  m -)  =  2 m2  +1,  m3  =  4 + 1  , 

^  =  6/77? +21/7^ +3  (3.11) 

As  the  closed  form  solution  covers  only  limited  cases,  it  has  been  shown  that  the 
most  efficient  method  of  determining  the  dynamic  range  is  by  using  a  search 
algorithm  [12,  13]. 
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B.  THE  RSNS  INTERFEROMETER  DESIGN 


The  design  of  an  RSNS  interferometer  has  been  simplified  into  eight  steps 
by  Lee  in  his  Master’s  thesis  [3].  The  eight  steps  are: 


1.  Determine  N,  the  number  of  channels  required. 

2.  Identify  an  integer  valued  modulus  for  each  channel  (mi ... mN ). 

3.  Determine  the  system  dynamic  range,  M ,  based  on  the 
chosen  moduli. 

4.  Determine  the  number  of  folds,  nf- ,  for  each  modulus. 


5. 


6. 


The  number  of  folds  is  determined  by 

M  M 


ni  = 


(3.12) 


PRSNS  2Nmi 

Determine  the  required  antenna  spacing,  dj ,  for  each  channel. 


.  am; i 

di  =  ni~  = 


2  4Nm, 

Re-mapping  the  Field  of  View  (FOV). 


(3.13) 


This  is  an  optional  step.  The  gain  of  an  antenna  reduces  as  the 
AOA  increases  from  the  broadside  and  to  get  a  high  signal-to- 
noise  ratio  (SNR)  for  digital  processing,  the  FOV  can  be 
compressed  into  a  suitable  range.  This  increases  resolution  and 
decreases  the  antenna  spacing  required  [4],  By  keeping  k  and  Vout 
constant  in  Equation  (2.5),  we  get 

Qfsin(6)max)  =  d  sin(6'max)  (3.14) 

where  <9max  is  the  maximum  FOV  and  d  is  the  antenna  spacing  of 
the  original  configuration.  The  quantities  6^ax  and  d'  are  the 

maximum  FOV  and  antenna  spacing  for  the  scaled  configuration. 
The  scale  factor,  is  defined  as 
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(3.15) 


g  _  ^  _  sin(^max ) 
d  s'n(^max) 

Combining  Equation  (3.13)  and  Equation  (3.15),  the  scaled  antenna 
spacing  can  be  determined. 

d'  =d/£  =  -^— £  (3.16) 

max  ,h  4/77,-A/ h  v  ' 

7.  Determine  the  thresholds. 


For  a  digital  architecture,  the  thresholds  to  amplitude  analyze  the 
folding  phase  waveform  into  the  RSNS.  It  has  been  shown  [4]  that 
for  a  normalized  folding  waveform,  the  threshold  for 
the  y'th  comparator  of  modulus  channel  rrij  can  be  determine  by 


vj,mi  =cos 


.  1 
m;  -  I  +  - 
'  7  2 

m. 


TV 


1  <  j  <mj. 


(3.17) 


V  J 

When  the  waveform  exceeds  a  threshold,  it  will  be  encoded  to  the 
next  integer  value.  Figure  7  shows  how  the  folding  waveform  in 
Table  1  can  be  encoded  by  Equation  (3.17).  At  the  boxed  in 


section  of  M,  the  symmetrical  residues  are  s1=0(/t?1=3)  and 


s2  =  3  (m2  =  4). 
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Figure  7.  Mapping  RSNS  folding  waveforms  to  threshold  values  (From  [5]). 


8.  Calculate  the  phase  adjustment  term,  £ ,  for  each  channel. 


From  Figure  7,  both  moduli  do  not  have  symmetry  about  the  center 
of  the  dynamic  range  at  n= 7.  A  phase  adjustment  term,^ ,  is  thus 
required  to  shift  the  center  of  the  dynamic  range  to  the  signal  AOA 
at  the  broadside.  From  Equation  (2.5),  the  phase  adjustment  term 
is  added: 

Vout(0)  =  cos(kds\n0  +  £)  (3.18) 

The  relationship  between  the  RSNS  coding  and  the  AOA  has  been 
obtained  at  this  point.  Using  this  relationship,  a  transfer  function  of 
the  estimated  AOA  against  the  true  AOA  can  be  plotted.  Previous 
studies  have  shown  that  the  estimated  AOA,  0'u ,  for  the  uth  bin  is 


given  as  [4]: 


d'  =  sin  1 


2u  +  1 


and  the  resolution,  ru ,  for  the  uth  bin  is  given  by: 


ru  =  sin 


-1 


2u-M  +  2 


-sin 


-1 


2n-M_ 


(3.19) 


(3.20) 
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c. 


DESIGN  EXAMPLE 


A  three-channel  RSNS  DF  system  was  designed  and  implemented  by 
Lee  in  his  master’s  thesis  [3].  Using  the  eight  steps  as  previously  described,  his 
design  is  used  as  an  example  to  explain  how  these  steps  are  carried  out  in  detail. 

1.  Determine  N,  the  number  of  channels  required. 

It  was  decided  a  three  channel  DF  system  was  to  be  designed.  Hence, 
N=3  and  a  total  of  four  antenna  elements  would  be  required. 

2.  Identify  an  integer  valued  modulus  for  each  channel 

The  moduli  set  of  [7  15  29]  was  chosen  for  high  resolution  as  well  as  a 
relatively  short  antenna  baseline.  It  is  important  to  note  that  there  will  be 
other  viable  alternatives  to  the  moduli  set  and  the  current  set  was  chosen 
because  it  consists  of  pairwise  relatively  prime  (PRP)  numbers. 

3.  Determine  the  system  dynamic  range,  M,  based  on  the 
chosen  moduli. 

Using  a  MATLAB  search  program,  the  dynamic  range  for  the  moduli  set 
was  found  to  b eM  =  444. 


4.  Determine  the  number  of  folds,  nf- ,  for  each  modulus. 

Applying  Equation  (3.12),  the  number  of  folds  for  each  channel  sequence 
was  found  to  be, 


nl  = 


n2  = 


n3  = 


M  _  444  _ 
2Nm<\  ~  2(3 )(7)  “ 

M  _  444 

2Nm2  ~  2(3)(15) 

M  _  444 

2/Vm-i  -  2(3)(29) 


10.57 


=  4.93 


=  2.55 
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5.  Determine  the  required  antenna  spacing,  dj ,  for  each  channel. 

At  this  step,  it  was  decided  that  a  signal  frequency  of  2.4  GHz  be  selected 
as  there  were  readily  available  component  hardware  in  the  laboratory, 
which  could  be  used  for  building  the  array.  For  a  signal  of  2.4  GHz,  the 
wavelength,  X,  is  12.5  cm.  Applying  Equation  (3.13),  the  distances 
between  pairs  of  antenna  elements  are, 

7  712  57 

d-\=  n-\  —  =  10.57  —  =66.07  cm 
2  2  j 

X  ( 1 2  5  A 

d2  =  n2  -  =  4.93  —  =30.83  cm 
2  2  J 

7  712  57 

do  =  no  -  =  4.93  —  =15.95  cm. 

O  O  r\  r\ 

V  ^  J 

6.  Re-mapping  the  Field  of  View  (FOV). 

The  FOV  was  not  scaled  and  hence  the  scaled  factor  £=1  (±90°) . 

7.  Determine  the  thresholds. 

Each  channel  would  have  a  number  of  thresholds  corresponding  to  the 
chosen  modulus  number.  For  a  modulus  of  7,  there  would  be  7 
thresholds.  Similarly,  modulus  of  15  and  29  would  have  15  and  29 
thresholds  respectively.  The  threshold  levels  for  the  chosen  moduli  [7  15 
29]  are  shown  graphically  in  Figure  8  below  using  Equation  (3.17). 
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Figure  8.  Threshold  levels  for  moduli  set  [7  15  29]. 

8.  Calculate  the  phase  adjustment  term,  £ ,  for  each  channel. 

The  phase  adjustment  terms,  ,  were  similarly  found  by  a  MATLAB 
search  program.  The  phase  adjustment  terms  are  shown  in  Table  2. 


Phase  Adjustment  Term 

Channel  1,^ 

-1.4960  radians 

Channel  2,^2 

0.0698  radians 

Channel  3,^3 

0.8666  radians 

Table  2.  Phase  adjustment  term  for  moduli  set  [7  15  29]  (From  [3]). 
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D.  SIMULATION  RESULTS 

Using  Equation  (3.18),  the  folding  waveforms  for  all  the  channels  were  plotted  in 
Figure  9.  The  corresponding  encoded  waveforms  were  also  plotted  in  Figure  10. 
By  using  Equation  (3.19),  the  estimated  AOA  can  be  plotted  against  the  true 
AOA  as  shown  in  Figure  1 1 . 


Angle  of  Arrival  (degrees) 

Figure  9.  Folding  waveforms  with  phase  adjustment. 


Angle  of  Arrival  (degrees) 


Figure  10.  Phase  adjusted  folding  waveforms  after  encoding. 
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T ransfer  Function  for  Moduli  [7  15  29] 


Figure  1 1 .  Transfer  function  for  demonstration  array  using 
moduli  [7  15  29]  (From  [3]). 

E.  EXPERIMENTAL  RESULTS 

A  RSNS-based  DF  array  was  built  and  tested  using  same  design  as 
described  in  Section  D.  The  experimental  results  are  shown  in  Figure  12.  The 
experimental  results  showed  that  although  there  are  errors  appearing 
intermittently  throughout  the  AOA,  the  basic  feature  of  the  transfer  function 
remains  similar  to  the  ideal  transfer  function  as  shown  in  Figure  11.  The  errors 
are  likely  due  to  antenna  spacing  error  and  phase  errors  from  the  demodulator 
boards  [3]. 
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Figure  12.  Experimental  results  for  demonstration  array 
using  moduli  [7  15  29]  (From  [3]). 
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IV.  RSNS  INTERFEROMETRY  WITH  VIRTUAL  SPACING 


Virtual  spacing  for  RSNS  interferomery  was  first  proposed  by  Chen  for  a 
two-channel  system  and  it  was  later  extended  to  an  A/-channel  system  by  Lee  [3, 
6].  This  chapter  introduced  how  quadrature  demodulation  works  followed  by  an 
explanation  of  how  a  two-channel  interferometer  works.  Finally,  a  three-channel 
virtual  spacing  DF  system  is  presented. 

A.  QUADRATURE  DEMODULATION 

The  bandpass  signal  can  be  expressed  by  the  in-phase  component  and 
the  qaudrature  component  as 

xband(t)  =  xi  cos (a>ct)-xQ  sin (coct)  (4.1 ) 

where  x ,  is  the  in-phase  component, 

xQ  is  the  quadrature  component,  and 


<dc  =  2n fc,  fc  is  the  carrier  frequency. 


demodulation 

Figure  13.  Quadrature  type  demodulation  of  two-channel  interferometer 

(From  [3]). 
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The  output  voltage  from  the  antenna  elements  as  shown  in  Figure  13  can 
be  expressed  as 

V(p„,t)  =  Re(A(p„)eJ(a’‘,>')  (4.2) 

where  n  =  1,2,3  is  the  antenna  element  index,  pn  is  the  position  at  antenna 


element  n,  t  is  the  time.  A(pn)  is  defined  as 

A(p„)  =  \A(p„pmp"'l=Al(Pn)  +  jAQ(Pn)  (4.3) 

The  output  voltage  can  be  mixed  with  a  local  oscillator  signal  to  get  the  in-phase 
and  quadrature  terms 

V/iPrvt)  =  V(pn,t)cos(coct) 

1  1  (4.4) 

=  2  A(Pn )  +  2  tA(Pn  )cos(2  coct)-AQ(pn  )sin(2  coct)] 

VQ(Pn,t)  =  -V(Pn,t)s\n((oct) 

1  1  (4-5) 

=  --[Ai(pn)s\n(2act)  +  ^(pn  )cos(2®cf )]  +  -  A?(Pn ) 

After  low-pass  filtering,  the  baseband  terms  that  remain  are 

Vl(Pn  )  =  ^Al(Pn)  =  \ \A(Pn  )|cos(./.(pn ))  (4.6) 

VQ(pn )  =  \ Aq(P„  )  =  ^A(Pn  )|  Sin  («>(p„ ))  (4.7) 

For  a  plane  wave,  let  |A(pn)|  =  1.  After  normalizing  Equations  (4.6)  and 

(4.7),  the  in-phase  and  quadrature  terms  of  the  received  signals  by  antenna 
element  n  are 


Vln=V,(pn)  =  cos(<fi(pn))  (4.8) 

VQn  =  VQ(Pn )  =  ^H<P(Pn ))  (4-9) 


The  phase  differential  term  between  antenna  element  two  and  three  can 
be  found  by 


Afe  =  <f>(P2)  ~  0(P3) 


=  arctan 


Vb2 

V  VI2  ) 


\  fy  \ 


-arctan 


'Q3 

V  VI3  J 


=  /cc/2  sin  ^ 


(4.10) 


24 


Similarly,  the  phase  differential  term  can  be  found  for  any  other  antenna 
pair.  The  in-phase  and  quadrature  outputs  for  any  channel  using  antenna 
element  r  as  the  reference,  can  be  defined  as 

VfN{0)  =  cos  (kdN  sind)  =  cos(A^m ) 

=  cos  (0(pn  ))cos  (</>(pr ))  +  sin(^(pn  ))sin  (</>(pr ))  (4.11) 

=  VjnVjr  +  VQnVQr 

Vqn(&)  =  sin  (kdN  sind)  =  sin(A^m ) 

=  sin(^(pn  ))cos (<t>(pr ))  -  cos ))sin(^(pf ))  (4.1 2) 

=  VinV/r ~VQnVQr 

The  phase  adjustment  term,  ,  is  required  for  channel  /  so  as  to  map  the  center 

of  the  dynamic  range  to  the  broadside.  From  Equation  (3.18),  the  voltage  from 
channel  /  can  be  expressed  as 

V-(<9)  =  cos(/cd/sin((9)  +  4';)  (4.13) 

Rearranging  Equation  (4.13), 

VN(d)  =  cos(/fd/sin(6)))cos(^/)-sin(/rc//sin(6)))sin(^/)  (4.14) 

V,  (d)  =  V,i  (d)cos {Ci)-VQi  (d)sin (Ci)  (4.15) 

B.  TWO-CHANNEL  VIRTUAL  SPACING  INTERFEROMETRY 

For  a  two-channel  RSNS  interferometer  as  shown  in  Figure  13,  the  moduli 
pair  m- \  and  m2  is  chosen.  The  dynamic  range  of  the  pair  is  M .  By  combining  the 
terms  as  given  in  Equation  (4.12)  and  Equation  (4.13),  the  in-phase  and 


qradrature  term  of  the  virtual  spacing  can  be  obtained: 

V/l  =  VnV/2  -VqiVq2  =  cos(/c(d1  +  d2)sind)  =  cos(kd\  sind)  (4.16) 
Vqi  =  Vq-)V/2  +  V/-)\/q2  =  sin(/f(d1 +d2)sind)  =  sin(/rd:jsind)  (4.17) 
Vj2  =  VnV/2  +Vq^VQ2  =  cos(/c(d1  -d2)  sind)  =  cos(kd2  sind)  (4.18) 
Vq2  =Vq'\Vi2-Vi'\Vq2  =sin(/f(d-| -d2)sind)  =  sin(/(d2sind)  (4.19) 

Now  suppose  that  the  chosen  moduli  set  [m-\  m2]  is  given  by 

=  hm2  (4.20) 
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where  h  is  the  ratio  of  modulus  m- \  and  m2 . 

From  Equation  (4.13)  and  Equation  (4.14),  using  Equation  (3.16), 

AM 


d\  =  d<\  +  d2  = 
d2  =  d- |  —  d2 


4N 

AM 


4N 


f  1  1  ^ 

e  AM 

f  1  +  A?^ 

4  =  -TTT 

Vmi  m2y 

AN 

1  m1  J 

f  .  .  \ 

1  1 

,  AM 

%  =  —77 

vmi  m2j 

AN 

^  m1  J 

„  AM'  e, 
,  AM'  c, 

4  =  Jfw4 


(4.21) 

(4.22) 


This  gives  the  relationship  between  a  new  moduli  set  [m\  m'2]  with 

dynamic  range  M'  and  scale  factor <f.  Combining  Equation  (4.18)  and  Equation 
(4.19), 


1  +  /7 

1 \-h 


_  m2 

m\ 


-»  m’2=m\ 


1  +  h 


(4.23) 


A  suitable  virtual  moduli  set  [m\  m'2]  can  be  found  by  using  Equation  (4.20)  using 
the  ratio  h.  A  simple  example  illustrates  how  this  virtual  moduli  set  would  work. 

Consider  a  case  of  a  two-channel  RSNS  interferometer  using  moduli  set 
[8  16],  the  dynamic  range  for  this  set  would  be  M=  35  and  the  scale  factor  £  =  1 . 
Let  the  frequency  be  2.4  GHz  (A  =  12.5  cm ).  The  ratio, 

h  =  ^l  =  - §_  =  1 
m2  16  2 

The  antenna  spacings  calculated  from  Equation  (3.13)  are 


d1  = 
d2  = 


MA  35x12.5 


ANrri'i 

MA 


4x2x8 

35x12.5 


=  6.8cm 

=  3.4cm 


(4.24) 

(4.25) 


4Nm2  4x2x16 

The  resolution  at  broadside  is  given  by  Equation  (3.20).  For  AOA  =  0°  at 
broadside,  the  bin  number,  u,  corresponds  to  90. 
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.  _if 2u-M  +  2) 
ru= 90=sin  - - 

l  M4  J 

=  sin-1  f  2(902-35  +  2 
l  35(1) 


.  _if 2u-M^\ 

-sin  — ^ — 

l  Mt  J 

2^1  .  _if 2(90)- 35 

—  -sin  — - — - - 

J  {  35(1) 


(4.26) 


=  3.3C 


From  Equation  (4.20), 


f  1  +  h  ^  M  +  o-  j 

m2=mi(j7^J  =  /77i  =3mi 

V  2  J 


(4.27) 


There  could  be  several  modular  pairs  which  satisfy  Equation  (4.21)  and  Table  3 
shows  some  of  the  possible  pairs  and  their  corresponding  dynamic  ranges. 


m\ 

7 

8 

9 

10 

m'2 

21 

24 

27 

30 

M' 

45 

51 

57 

63 

Table  3.  Moduli  sets  and  dynamic  range  for  h  =  %  (From  [6]). 


Suppose  that  the  chosen  virtual  moduli  set  is  [m\  m'2]  =  [10  30]  with 

M'  =  63.  The  virtual  scale  factor  can  be  found  from  Equation  (4.18)  or  Equation 
(4.19), 


d\=d^+d2  =  (6.8  +  3.4)  =  1Z5(63)<f 
1  1  ^  4(2)(10) 


4'  =  1 .04 

The  resolution  of  the  virtual  moduli  set  is 


,  .  _if  2u-M'  +  2 )  ,_if  2 u-M' 

ru- 90=s|n  = - sin  — ^ 

1  M'4'  J  (  M'4' 


.  _if 2(90)-63  +  2  ]  .  _if 2(90) -63 

=  sin  — - — - -  -sin  — - — - - 

f  63(1.04)  J  {  63(1.04) 

=  1.8° 


(4.28) 


(4.29) 
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The  virtual  moduli  set  improved  upon  the  resolution  of  the  original  array  by 
almost  100%.  In  theory,  super  resolution  could  be  achieved  with  a  large  number 
virtual  moduli  set  [m\  m >2].  However,  the  larger  the  m\  or  m2 ,  there  would  be  a 

corresponding  increase  in  the  number  of  thresholds.  As  shown  in  Figure  8,  these 
levels  of  threshold  are  closer  together  at  the  peaks  and  troughs  of  the  folding 
waveform.  Having  the  threshold  levels  closer  can  lead  to  quantization  error  in 
the  presence  of  noise  if  the  signal-to-noise  ratio  is  low.  It  has  been  shown  by 
Chen  that  in  a  low  SNR  environment,  every  moduli  set  would  have  regions  where 
the  transfer  function  “fades,”  i.e.,  where  gross  errors  tend  to  occur  [6]. 

C.  THREE-CHANNEL  VIRTUAL  SPACING  INTERFEROMETRY 


The  virtual  RSNS  DF  can  be  similarly  extended  to  a  A/-channel  system  [3]. 
Let  us  consider  a  case  of  three-channel  virtual  DF  array  where  the  moduli  set  are 
related  by 

/T7-)  =  hm2  =  h2m2  (4.30) 

The  ratio  h  must  satisfy  the  condition  h  <  ^  for  three  integer  values  of  m,-  to 

exist.  Similar  to  the  two-channel  system,  the  virtual  spacing  may  be  extended  to 
a  three-channel  system  by  defining  d3  as 


do  =  d- i  —  do  —  do  = 


AM 


%  f 


AN 


J 1_ 

m1 


J _ 1_ 

m2  m3 


4  = 


AM 


AN 


1  -h-h‘ 
m1 


AM' 

ANmy  (4.31) 


Rearranging  Equation  (4.18)  with  Equation  (4.19)  and  Equation  (4.30),  we  get 

+  h^ 


1  -h 


m2 

m\ 


(4.32) 


^  1  +  h  ^  _m2 

vl -h-h2)  m\ 


(4.33) 


Depending  on  the  ratio  h,  it  is  possible  to  select  integer  value  of  [m\  m2  m2]  so 

as  to  form  other  virtual  moduli  sets.  A  simple  algorithm  can  be  implemented  to 
post-process  the  AOA  information  using  the  virtual  moduli  set  to  increase  the 
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resolution  of  the  DF  array.  A  simple  algorithm  flow  with  two  levels  of  virtual 
resolution  processing  is  shown  in  Figure  14. 


•  Estimate  AOA  (eAOAl)  from  moduli  set 

•  Estimate  virtual  AOA  (eAOA2  and  eAOA3)  for 
2  set  of  virtual  moduli  where  resolution  of 
eAOA2>resolution  of  eAOA3 


Figure  14.  Algorithm  for  two  levels  of  virtual  resolution  processing  (After  [3]). 


29 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


30 


V.  DESIGN  CONSIDERATIONS  FOR  THREE-CHANNEL  RSNS 
VIRTUAL  SPACING  DIRECTION  FINDING  SYSTEM 


This  chapter  presents  some  of  the  considerations  in  designing  a  three- 
channel  RSNS  virtual  spacing  DF  system. 

A.  LIMITED  REAL  MODULI  SET 

Let  [m-i  m2  m3]  be  defined  as  the  real  moduli  set  where  m1  <  m2  <  at?3 .  In 

order  for  virtual  resolution  to  work,  the  real  moduli  set  must  satisfy  Equation 
(4.27).  Let  the  ratio,  R,  be  defined  as 

r=4  m 

h 2 

As  mentioned  in  the  previous  chapter,  h<\ ^  is  required  for  three  distinct 

positive  integer  values  of  the  moduli  set  to  exist.  From  Equation  (5.1),  R>  2 
follows. 

Based  on  these  conditions,  some  of  the  possible  real  moduli  sets  which 
can  be  selected  for  a  three-channel  RSNS  DF  system  are  shown  in  Table  4 
below. 
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[m1  m2  m3] 

M 

h  =  y2,  R  =  2 

[2  4  8] 

29 

[3  6  12] 

41 

[4  8  16] 

53 

[5  10  20] 

65 

h  =  y2,  R  =  3 

[2  6  18] 

59 

[3  9  27] 

86 

[4  12  36] 

113 

[5  15  45] 

140 

h=y4,R=4 

[2  8  32] 

101 

[3  12  48] 

149 

[4  16  64] 

197 

[5  20  80] 

245 

Table  4.  Real  moduli  set  [m1  m2  m3]  and  dynamic  range  for 

several  values  of  h. 


B.  LIMITED  VIRTUAL  MODULI  SET 


Let  [m\  m'2  m3]  be  defined  as  the  virtual  moduli  set  where  m\  <  m'2  <  m3 . 
For  the  virtual  moduli  set  to  have  positive  integer  value,  it  needs  to  satisfy 
Equation  (4.32)  and  Equation  (4.33).  Let  the  ratios,  R-\  and  R2,  be  defined  as 


R<  = 


m2 


Ro  = 


m\ 

m'3  f 


m\ 


1  +  h 
1  -h 
1  +  h 


U-h-h' 


(5.2) 

(5.3) 


For  the  virtual  moduli  set  to  have  positive  integer  value,  m\  must  satisfy  both 
Equation  (5.2)  and  Equation  (5.3).  Hence  m\  has  to  be  the  multiple  of  lowest 

common  denominator  of  R-\  and  R2  ■  For  h  =  ^ ,  R^=  3  and  R2  =  6 .  Since  the 
lowest  common  denominator  of  R 1  and  R2  is  1 ,  there  is  no  restriction  on  the 
value  of  m\.  For  h  =  ^,  R1=2  and  R2=^%-  As  the  lowest  common 
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denominator  of  R- |  and  R2  is  5,  m\  can  be  any  multiple  of  5.  Table  5  lists  some 
of  possible  the  virtual  moduli  sets  for  selected  ratios. 


[m\  m2  m3] 

M' 

h~y2 

[2  6  12] 

41 

[3  9  18] 

59 

R|  =  2,  R2  =  6 

[4  12  24] 

77 

[9  27  54] 

167 

h=Y 

[5  10  12] 

89 

[10  20  24] 

179 

R,=2,  R2='i% 

[15  30  36] 

269 

[20  40  48] 

359 

h=yA 

[33  55  60] 

853 

[66  110  120] 

1708 

[99  165  180] 

2563 

[132  220  240] 

3418 

Table  5.  Virtual  moduli  sets  [m\  m'2  m3  ]  and  dynamic  range  for 

several  values  of  h. 


C.  SENSITIVITY  OF  SIZE  OF  ANTENNA 

The  size  of  an  antenna  array  is  determined  by  Equation  (3.16).  If  the 
design  frequency  f  is  fixed  at  2.4  GHz,  the  antenna  size  is  only  dependent  on 

the  selection  of  modulus,  m,-,  dynamic  range  M  and  scale  factor  Since  the 
dynamic  range  M  is  determined  by  modulus  m;- ,  the  antenna  size  is  only 
dependent  on  two  parameters,  m,-  and  £  ■ 

The  antenna  spacing  of  the  moduli  set  is  listed  in  Table  6.  It  shows  that 
the  antenna  array  spacing  is  relatively  insensitive  to  the  selection  of  the  real 
moduli  set  and  more  dependent  on  the  ratio  h. 
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[m1  m2  m3] 

d- |  (cm) 

d2  (cm) 

d3  (cm) 

h~y2 

[2  4  8] 

15.1 

7.6 

3.8 

[3  6  12] 

14.2 

7.1 

3.6 

4  =  1 

[4  8  16] 

13.8 

6.9 

3.5 

[5  10  20] 

13.5 

6.8 

3.4 

h=X 

[2  6  18] 

30.7 

10.2 

3.4 

[3  9  27] 

29.9 

10.0 

3.3 

4=1 

[4  12  36] 

29.4 

9.8 

3.3 

[5  15  45] 

29.2 

9.7 

3.2 

h-/4 

[2  8  32] 

52.6 

13.2 

3.3 

[3  12  48] 

51.7 

12.9 

3.2 

£  =  1 

[4  16  64] 

51.3 

12.8 

3.2 

[5  20  80] 

51.0 

12.8 

3.2 

Table  6.  Antenna  spacing  for  [m- \  m2  m3] 
(scale  factor  £  =  1  and  frequency  at  2.4  GHz). 


Changing  the  scale  factor  to  £  =  1.154  reduces  the  FOV  to  ±60°  but  it 
increases  the  array  spacing  slightly  as  shown  in  Table  7. 


[m-|  m2  m3] 

d-|  (cm) 

d2  (cm) 

d3  (cm) 

h=y2 

[2  4  8] 

17.4 

8.7 

4.4 

[3  6  12] 

16.4 

8.2 

4.1 

£  =  1.154 

[4  8  16] 

15.9 

8.0 

4.0 

[5  10  20] 

15.6 

7.8 

3.9 

h  =  Y$ 

[2  6  18] 

35.5 

11.8 

3.9 

[3  9  27] 

34.5 

11.5 

3.8 

£  =  1.154 

[4  12  36] 

34.0 

11.3 

3.8 

[5  15  45] 

33.7 

11.2 

3.7 

"-y* 

[2  8  32] 

60.7 

15.2 

3.8 

[3  12  48] 

59.7 

14.9 

3.7 

£  =  1.154 

[4  16  64] 

59.2 

14.8 

3.7 

[5  20  80] 

58.9 

14.7 

3.7 

Table  7.  Antenna  spacing  for  [m- \  m2  m3] 
(scale  factor  £  =  1.154  and  frequency  at  2.4  GHz).. 
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D.  FIELD  OF  VIEW  LIMITATION 

The  FOV  of  the  real  moduli  set  is  determined  by  the  Equation  (3.15)  while 
the  FOV  of  the  virtual  moduli  is  determined  by  Equation  (4.21),  Equation  (4.22)  or 
Equation  (4.31).  The  critical  parameter  in  determining  the  FOV  is  the  ratio  h. 
Table  8  and  Table  9  list  the  FOVs  of  selected  combinations  of  real  and  virtual 
moduli  sets  for  £  =  1  and  £  =  1.154. 


[m1  m2  m3] 

[m\  m'2  m'3] 

FOV(deg) 

FOV'(deg) 

\CN 

II  II 

■C  HP 

[3  6  12] 

[4  12  24] 

±90 

±69.9 

[9  27  54] 

±90 

±64.8 

[4  8  16] 

[4  12  24] 

±90 

±75.6 

[9  27  54] 

±90 

±69.0 

[5  10  20] 

[4  12  24] 

±90 

±80.8 

[9  27  54] 

±90 

±72.1 

\00 

II  II 

[2  6  18] 

[5  10  12] 

±90 

±26.9 

[10  20  24] 

±90 

±27.1 

[3  9  27] 

[5  10  12] 

±90 

±27.8 

[10  20  24] 

±90 

±27.9 

[4  12  36] 

[5  10  12] 

±90 

±28.2 

[10  20  24] 

±90 

±28.4 

II  II 

■C  Hn 

[2  8  32] 

[33  55  60] 

±90 

±24.2 

[66  110  120] 

±90 

±24.2 

[3  12  48] 

[33  55  60] 

±90 

±24.6 

[66  110  120] 

±90 

±24.6 

[4  16  64] 

[33  55  60] 

±90 

±24.8 

[66  110  120] 

±90 

±24.9 

Table  8.  Field  of  view  for  selected  moduli  set,  £  =  1 . 
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[m\  m2  m3] 

[m\  m'2  m'3] 

FOV(deg) 

FOV'(deg) 

h~Yi 

£  =  1.154 

[3  6  12] 

[4  12  24] 

±60 

±54.4 

[9  27  54] 

±60 

±51.6 

[4  8  16] 

[4  12  24] 

±60 

±57.0 

[9  27  54] 

±60 

±54.0 

[5  10  20] 

[4  12  24] 

±60 

±58.8 

[9  27  54] 

±60 

±55.5 

h  =  Y$ 

£  =  1.154 

[2  6  18] 

[5  10  12] 

±60 

±23.1 

[10  20  24] 

±60 

±23.2 

[3  9  27] 

[5  10  12] 

±60 

±23.8 

[10  20  24] 

±60 

±23.9 

[4  12  36] 

[5  10  12] 

±60 

±24.2 

[10  20  24] 

±60 

±24.3 

h=y4 

£  =  1.154 

[2  8  32] 

[33  55  60] 

±60 

±20.8 

[66  110  120] 

±60 

±20.8 

[3  12  48] 

[33  55  60] 

±60 

±21.1 

[66  110  120] 

±60 

±21.2 

[4  16  64] 

[33  55  60] 

±60 

±21.3 

[66  110  120] 

±60 

±21.3 

Table  9.  Field  of  view  for  selected  moduli  set,  £  =  1 .1 54 . 

The  DF  system’s  FOV  is  determined  primarily  by  h  and  is  relatively 
insensitive  to  the  selection  of  the  moduli  set.  As  the  ratio  h  decreases,  the 
system’s  FOV  reduces  accordingly.  For  a  given  operational  requirement  of  an 
instantaneous  FOV,  the  choice  of  the  ratio  h  could  be  limited. 

E.  SELECTION  OF  MODULUS  RATIO 

Let  the  phase  differences  between  every  combination  of  element  of  a 
three-channel  interferometer  be  defined  as  indicated  in  Figure  15.  An  ambiguity 
diagram  can  be  used  to  determine  which  h  ratio  should  be  chosen.  An  ambiguity 
diagram  is  a  phase-plane  plot  comparing  the  phase  difference  of  an 
interferometer  [14].  The  closer  the  lines  on  the  ambiguity  diagrams,  the  poorer 
the  interferometer  will  perform  in  determining  the  actual  phase  difference.  To 
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further  explain  how  the  ambiguity  diagram  works,  consider  a  simple  three 
element  interferometer  with  d1  =  and  d2  =  2^j/. ^ j .  The  phase-plane  plot  of 

the  corresponding  ^and  ^2is  given  in  Figure  16. 


4  3  2  1 

V  V  V  V 


* — 1 

_ 4 _ . 

- 4 


Figure  15.  Phase  difference  for  a  three-channel  interferometer. 


Figure  16.  Phase-plane  plot  for  d,  =3(^ )  and  d2  =2(/^)  (From  [14]). 
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Since  ^and  (j>2  can  be  known  only  to  within  2 n  radian,  Figure  16  can  be 
normalized  to  as  shown  in  Figure  17.  The  ambiguity  diagrams  for  selected 
ratios  of  /?  =  ^ ,  h  =  ^  and  h  =  ^  are  given  in  Figures  18  to  20.  These  figures 

showed  that  the  higher  h  ratio  gives  the  best  performance  in  resolving  ambiguity 
between  the  phases  as  the  lines  are  further  apart. 


Figure  17.  Ambiguity  diagram  for  d1  =3(^2)  and  d2  =  2^^j  (From  [14]). 

In  designing  a  RSNS  virtual  spacing  DF  system,  it  is  important  to  select 
the  highest  possible  h  ratio.  For  any  given  h  ratio  and  frequency,  the  size  of  the 
antenna  array  as  well  as  the  FOV  of  the  system  would  not  change  significantly 
once  the  ratio  has  been  selected  The  highest  possible  h  ratio  gives  the  best 
performance  in  resolving  phase  ambiguity. 
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Figure  1 8.  Ambiguity  diagram  for  h  =  ^ 
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Figure  1 9.  Ambiguity  diagram  for  h  =  ^ . 
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Figure  20.  Ambiguity  diagram  for  h  =  ^ . 
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VI.  MODELING  AND  SIMULATION 


To  understand  and  predict  the  performance  of  a  three-channel  RSNS 
virtual  spacing  DF  system,  simulations  in  MATLAB  were  carried  out.  This 
chapter  introduces  all  the  models  that  were  used  in  the  MATLAB  simulations. 

A.  NOISE  MODEL 

It  is  assumed  that  the  noise  in  the  system  is  random  in  nature  and 
Gaussian  distributed.  The  bandpass  Gaussian  random  noise  signal  can  be 
represented  by 

W(t)  =  Re(jA{t)\ej<p{t)ejc°ct]j  (6.1) 

where 

|A(f)|  is  random  process  with  Rayleigh  first-order  probability  density 
function;  and 

q o(t)-  random  process  with  uniform  first-order  probability  over  range 

(0,2*). 

Using  Euler’s  formula,  Equation  (6.1)  can  be  written  as: 

W(t)  =  W,  cos(coct)-WQ(t)s\n(coct)  (6.2) 

where 

lV,(()  =  |A(t)|cos(^())  (6.3) 

Wb<f>  =  Msin(p(f))  (6.4) 

From  these  equations  and  combining  Equation  (4.4)  and  Equation  (4.5), 
the  Gaussian  noise  model  can  be  represented  by 

\ fehe(Pn.0  =  V(p#J,f)  +  IV/(0  (6.5) 

VQNniJPrvt)  =  V(Pn,t)  +  WQ(t)  (6.6) 
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B.  MONTE  CARLO  SIMULATION 


Monte  Carlo  methods  are  commonly  used  in  simulating  physical  or 
mathematical  systems  when  a  closed  form  statistical  solution  is  too  complex  to 
resolve.  It  uses  the  Monte  Carlo  algorithm  which  can  be  easily  implemented  with 
a  computer  to  model  random  processes.  The  Monte  Carlo  simulation  is  carried 
out  to  determine  the  performance  of  the  three-channel  RSNS  virtual  spacing  DF 
system  with  random  white  Gaussian  noise.  First,  a  review  of  the  Monte  Carlo 
algorithm  is  given  below  [15]. 


Let  a  random  variable  X  consist  of  both  the  signal  voltage  V  together  with 
random  noise  W,  centered  on  zero  with  variance  of  one. 


X  =  V  +  W 

The  probability  of  X>z  for  any  given  V  where  z  is  constant  is 


P(z)  =  P 


f  z  ^ 
X>  — 
v  Vj 


(6.7) 

(6.8) 


By  generating  a  series  of  statistically  independent  and  identically 
distributed  Gaussian  random  variable  w(n),  where  n=1,2,...,A/  times,  Equation 
(6.8)  can  be  estimated.  A  new  set  of  random  variables,  x(n),  can  be  found  by 
adding  a  constant  vto  the  w(n). 


x(n)  =  v  +  w(n) 

Let  a  new  random  variable  y(n)  be  defined  as 

fO,  if  x(n)<z 


y(n)  = 


1,  if  x(n)>z 


(6.9) 

(6.10) 


The  estimation  of  Equation  (6.8)  is 

1  N 

P{z)  =  -Yy{n) 


(6.11) 


The  mean,  E(*},  and  variance,  cr2 


E(P(z)) 


,  of  P(z)  is 

1  N  1 
=jji:E(y("))  =  HNpM=pM 


(6.12) 
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CT/S<z)  =  E(p<z>-  p<z>)  =  e(p2(z)}-  P2(Z)  =  lp(z)(1  -P(z)) 
Let  the  required  outcome  probability  of  the  estimation,  a ,  be  defined  by 

a<p(|P(z)-P(z)|<s) 

where  the  relative  error  s  is  defined  as, 

P(z)-P(z) 

P(z) 


Using  Chebyshev’s  inequality, 


p(|P(z)-P(z)  <  >  1 


Using  Equation  (6.13), 


1 


N  > 


P(z)(1-P(z» 
(1  -a)z2 


(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 


Equation  (6.17)  determines  the  number  of  times  the  simulation  must  be  run  to 
satisfy  Equation  (6.14)  for  any  given  a  and^.  The  Monte  Carlo  algorithm  is 
summarized  by  the  flowchart  in  Figure  21 . 


45 


Figure  21 .  Monte  Carlo  algorithm. 

C.  HARDWARE  SETUP 

The  simulations  were  based  on  the  hardware  used  by  Lee  in  his  three- 
channel  DF  antenna  array  [3].  Figure  22  shows  the  block  diagram  of  all  the 
hardware  used  in  the  array. 
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Ground  Plane 


Figure  22.  Hardware  block  diagram  of  a  three-channel  DF  antenna  array. 


To  determine  the  SNR  of  the  I  and  Q  data  after  the  demodulator  board,  a 
calculation  of  the  noise  figure  was  done  for  one  receiver  channel.  The 
specifications  of  the  hardware  used  are  given  in  Figure  23. 
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Figure  23.  Specifications  for  antenna  array’s  noise  figure. 
The  receiver’s  noise  figure  as  shown  in  the  dashed  box  is  given  by 


Fe=F- ,+- 


1  Fo-1 

-  +  — —  +  - 


Fa  -1 


G-)G2  G-1G9G- 


(6.18) 
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where 


F-|  =  yG^  =  cable  loss  from  antenna  to  LNA  =  0.5  dB, 

F2  =  noise  figure  of  LNA  =  1.7  dB, 

G2  =gain  of  LNA  =  23  dB, 

F3  =  yG^  =  cable  loss  from  LNA  to  AD  8347  demodulator  =  0.65  dB, 
F4  =  noise  figure  of  AD  8347  demodulator  =11.7  dB, 


G4  =total  gain  of  the  demodulator  baseband  and  RF  amplifier. 


From  Equation  (6.18), 


_  „  1.48-1  1.16-1 

FQ  =1.12  + - + - 


14.79-1 


0.89  0.89x199.5  0.89x199.5x0.86 

Fe=  1-77 

Te  =(Fe-1)ra  =(1.77-1)290  =  223  K 
Ts=Te+Ta  =223  +  300  =  523  K 


The  noise  power  is 


(6.19) 

(6.20) 
(6.21) 
(6.22) 


N0  =  kBTsB 

=  1 .38  x  1 0-23  x  523  x  50  x  1 06  /e  oox 

(6-23) 

=  3.61x10-13 
=  -124.4  dB 

where  B  is  the  instantaneous  receiver  bandwidth  (50  MHz).  The  minimum  RF 
power  that  can  be  input  to  the  AD  8347  depends  on  the  sensitivity  of  the 
demodulator.  From  the  specification,  it  ranges  from  -65  dBm  to  0  dBm. 
Therefore  the  minimum  power  that  the  antenna  dipole  can  receive  is 

Prec  min  =  -65  +  0.65  -  23  +  0.5  -  6.55  =  -93.4  dBm  (6.24) 

The  minimum  SNR  output  from  the  AD8347  demodulator  board  is 

SA/Rmin  =-93.4 -(-124.4)  =  31.0  dB  (6.25) 

Assuming  that  atmospheric  absorption  is  negligible,  for  any  isotropic 
emitter  transmitting  an  average  power  of  Pt  at  distance  R,  the  SNR  can  be 
determined  by 
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SNR  = 


(6.26) 


Pt 

4ttR2N0 

Based  on  the  calculated  SNR,  the  variance  of  the  noise  is  determined  from 

v} 

SNR  =  ^r  (6.27) 

2cr2 

where  V/  is  the  channel  output  voltage  defined  in  Equation  (4.16).  For  a  one 
watt  emitter,  the  maximum  detection  range  of  the  current  array  is 

SNRmin= - ^4—  (6.28) 

1 259  = - ^ - txt  (6.29) 

4 ffxRma/  x  3.61x10 

Rmax  =13.2  km  (6.30) 


D.  ANTENNA  PATTERN 

The  antenna  array  used  by  Lee  utilizes  vertical  printed  circuit  dipoles  and 
was  designed  for  a  nominal  frequency  of  2.4  GHz.  It  has  been  tested  in  the 
anechoic  chamber  by  Burgstaller  to  determine  the  dipole  antenna  pattern  [16]. 
Figure  24  shows  the  measured  element  patterns  for  a  dipole  antenna  at  2.4  GHz. 
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Figure  24.  Dipole  element  pattern  measured  at  2.4  GHz  (From  [13]). 


The  gain  of  the  dipole  element  was  measured  to  be  6.55  dB  or  4.52  at 
broadside.  The  array  estimates  the  AOA  in  the  horizontal  plane  (H-plane).  For 
the  simulation,  a  simple  cosine  model  was  used  to  replicate  the  antenna  pattern 
for  various  AOA 

Ga  =  4.52xcos(#)  (6.31) 

All  the  tools  used  in  the  simulation  have  been  presented  and  the  next 
chapter  presents  the  results  from  simulations  that  were  carried  out. 
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VII.  RESULTS  AND  ANALYSIS 


With  all  the  tools  of  the  simulation  in  place,  a  sensitivity  analysis  was 
carried  out  to  find  out  how  some  of  these  errors  might  affect  the  performance  of 
the  RSNS-based  DF  system. 

A.  BASELINE  DESIGN  PARAMETER 

A  three-channel  RSNS  virtual  spacing  DF  design  was  chosen  as  a 
baseline  design  before  any  sensitivity  analysis  could  be  carried  out.  Using  the 
steps  described  in  Chapter  III,  the  design  parameters  for  a  three-channel  RSNS 
virtual  DF  system  were  chosen.  It  was  decided  that  the  design  would  be  based 
on  the  hardware  shown  in  Figure  22;  hence  the  design  frequency  of  2.4  GHz  was 
selected. 

1.  Determine  N,  the  number  of  channels  required. 

A  three-channel  DF  system  was  to  be  designed.  Hence,  A/=3  and  a 
total  of  four  antenna  elements  would  be  required. 

2.  Identify  an  integer  valued  modulus  for  each  channel  (mi,  m2 

m3). 

Before  identifying  the  moduli  sets,  the  ratio  h  must  first  be  selected. 
The  three-channel  RSNS  virtual  spacing  DF  system  in  Lee’s  thesis 
was  based  on  h  =  y2 .  The  size  of  the  antenna  array  built  was 

about  14  cm.  Experiment  results  shown  that  a  number  of  gross 
errors  appeared  throughout  the  transfer  function.  Although  the 
sources  of  errors  were  not  fully  investigated  by  Lee,  it  was 
suspected  that  the  errors  were  due  to  mutual  couplings  as  the 
antenna  spacing  between  the  antenna  elements  was  small.  It  was 
decided  that  h  =  %  would  be  appropriate  to  reduce  the  effects  from 

mutual  couplings  as  well  as  to  obtain  a  reasonable  system’s  FOV. 
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Based  on  h  =  % ,  the  real  moduli  set  of  [2  6  18]  was  selected  so  as 

to  reduce  quantization  error  when  noise  is  present.  The  virtual 
moduli  sets  of  [5  10  12]  and  [15  30  36]  were  selected  to  enable  the 
DF  system  to  have  sufficient  resolution. 

3.  Determine  the  system  dynamic  range,  M,  based  on  the 
chosen  moduli. 

From  Table  4  and  Table  5,  the  dynamic  ranges  for  the  all  the 
moduli  sets  are  listed  in  Table  10. 


Moduli  set 

Dynamic  Range 

[2  6  18] 

M Real  =  89 

[5  10  12] 

M Virtual!  =  89 

[15  30  36] 

^Virtual  2  =  269 

Table  10.  Dynamic  ranges  for  baseline  design. 


4.  Determine  the  numbers  of  folds,  nf- ,  for  each  modulus. 

The  numbers  of  folds  are  given  in  Table  11. 
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Moduli  Set 

Number  of  folds 

n1  =  4.92 

[2  6  18] 

n2  =  1 .64 

n3  =  0.55 

n\  =  2.97 

[5  10  12] 

n'2  =  1 .48 

n3  =  1 .24 

n\  =  2.99 

[15  30  36] 

ri2  =  1 .49 

ri3  =  1 .25 

Table  1 1 .  Numbers  of  folds  for  baseline  design. 


Determine  the  required  antenna  spacing,  dj ,  for  each  channel. 


The  required  antenna  spacing  for  all  moduli  sets  are  given  below  in 
Table  12. 


Moduli  Set 

Antenna  spacing 

[2  6  18] 

d- 1  =  30.73  cm 
d2  =10.24  cm 
d3  =  3.41  cm 

[5  10  12] 

d\  =  39.33  cm 
d2  =19.67  cm 
d3  =16.39  cm 

[15  30  36] 

d\  =  39.33  cm 
d2  =19.67  cm 
d3  =16.39  cm 

Table  12.  Antenna  spacing  for  baseline  design. 
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6. 


Re-mapping  the  Field  of  View  (FOV). 

The  real  moduli  set’s  FOV  was  not  scaled  and  hence  the  scale 
factor  %reai  =  1 .  The  FOV  of  the  three-channel  system  is  given  in 
Table  13. 


Table  13.  Scale  factor  and  FOV  for  baseline  design. 


7.  Determine  the  thresholds. 

The  threshold  levels  for  all  the  moduli  sets  are  shown  graphically  in 
Figure  25-27. 


/7T|  =  2  A7?2  =  6 


Figure  25.  Threshold  levels  for  real  moduli  set  [2  6  18]. 
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8. 


Calculate  the  phase  adjustment  term,  47,  for  each  channel. 


The  phase  adjustment  terms,  47  >  for  the  moduli  sets  are  listed  in 
Table  14. 


Moduli  set 

Phase  Adjustment  term,  C, 

47  =  2.36 

[2  6  18] 

4-2  =-1-48 

473=1.54 

4'j  =  -0.94 

[5  10  12] 

42  =-2.15 

a  =1-79 

4'j  =  -1 .15 

[15  30  36] 

42  =-2.18 

4-^  =1.82 

Table  14.  Phase  adjustment  term  for  baseline  design. 


The  ideal  transfer  function  for  the  baseline  design  is  given  in  Figure  28. 
The  baseline  design  has  a  broadside  resolution  of  0.19  deg  RMS.  In  practice,  it 
is  not  possible  to  achieve  this  accuracy  as  there  would  be  errors  which  have  not 
been  accounted  for  in  the  current  baseline  design.  A  one  watt  emitter  at  a 
detection  range  of  12  km  for  the  current  hardware  is  used  as  a  scenario.  This 
scenario  gives  a  SNR  of  approximately  38  dB.  Using  this  scenario  and  with  the 
addition  of  Gaussian  noise  as  described  in  Chapter  VI,  the  transfer  function  is 
given  in  Figure  29.  The  RMS  error  is  approximately  0.5  deg  RMS. 
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Figure  28.  Transfer  function  for  the  baseline  design. 


Figure  29.  Transfer  function  with  noise. 


The  transfer  function  shown  in  Figure  29  represents  a  snapshot  of  the 
transfer  function  in  time.  Over  time,  the  transfer  function  would  vary  due  to 
random  noise.  In  order  to  figure  out  how  the  noise  would  affect  the  probability  of 
estimating  the  correct  AOA  over  time,  a  Monte  Carlo  simulation  was  carried  out. 
For  the  Monte  Carlo  simulation,  the  following  parameters  were  used: 
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Relative  error,  s  =  3% . 


Required  outcome  probability  ,  a  =  90% . 

For  the  given  scenario,  there  is  a  greater  than  50%  probability  that  the 
correct  AOA  will  be  estimated  within  the  FOV  as  shown  in  Figure  30.  In  order  to 
satisfy  Equation  (6.14),  the  number  of  times  the  simulation  needs  to  be  run  is 
2777.  In  practice,  2777  readings  would  need  to  be  taken  from  the  emitter  when 
the  determining  the  AOA  for  the  results  of  the  Monte  Carlo  simulation  to  hold 
true.  As  this  is  not  realistc,  a  comparison  between  the  probabilities  of  taking  a 
reduced  sample  size  of  100  readings  against  2777  readings  is  given  in  Figure  31. 
It  shows  that  a  reduced  sample  of  100  readings  has  good  correlation  to  the  full 
Monte  Carlo  simulation  results.  Therefore,  even  if  a  reduced  sample  size  of  100 
readings  is  used,  the  results  from  the  Monte  Carlo  simulation  would  still  hold  true. 
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Figure  30.  Probability  of  correct  AOA  estimation  for 

baseline  design. 
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Figure  31 .  Probability  of  correct  AOA  estimation  for 
full  sample  size  vs  100  sample  size. 

B.  EFFECTS  OF  SPACING  ERROR 

As  it  is  unlikely  that  the  array  would  be  precisely  constructed,  the  effects  of 
spacing  error  were  investigated.  By  adding  a  manufacturing  error  of  2  mm  and 
4  mm  along  the  array  axis  into  the  element  spacing,  the  transfer  functions  were 
plotted  in  Figure  32.  The  spacing  error  does  not  affect  the  transfer  function 
significantly  for  spacing  errors  of  2  mm  while  spacing  errors  of  4  mm  introduce 
more  fading  regions  away  from  the  broadside.  These  spacing  errors  correspond 
to  0.016  and  0.032  wavelength  displacements  in  the  horizontal  direction.  Note 
that  position  errors  along  the  array  axis  do  not  affect  AOA  estimate  near  the 
broadside.  Position  errors  out  of  the  array  plane  will  affect  the  AOA  estimates 
near  the  broadside. 
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Figure  32.  Transfer  function  with  noise  and  spacing  errors. 

The  probability  of  correct  AOA  estimation  for  spacing  errors  of  2  mm  and 
4  mm  is  given  in  Figure  33.  Comparing  Figure  30  and  Figure  33,  spacing  errors 
do  not  affect  the  probability  of  correct  AOA  estimation.  These  results  suggest 
that  small  spacing  error  does  not  have  significant  effect  on  the  performance  of 
the  RSNS  virtual  spacing  DF  system. 


Figure  33.  Probability  of  correct  AOA  estimation  with  spacing  errors. 
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C.  EFFECTS  OF  PHASE  ERROR 

There  may  be  residual  phase  errors  introduced  into  the  DF  system  that 
may  not  be  possible  to  remove  during  the  calibration.  Such  phase  errors  could 
result  from  a  different  antenna  temperature  and  cable  phase  error.  By  adding  a 
phase  error  of  5  degrees  and  10  degrees  to  all  the  channels,  the  transfer  function 
is  plotted  in  Figure  34.  Comparing  Figure  29  with  Figure  34  shows  that  the 
phase  error  affected  the  transfer  function  accuracy  significantly.  As  the  phase 
error  increases,  the  transfer  function  has  more  fading  regions. 


Figure  34.  Transfer  function  with  noise  and  phase  errors. 

The  probability  of  correct  AOA  estimation  with  phase  errors  is  plotted  in 
Figure  35.  It  showed  with  phase  error  the  probability  of  correct  AOA  estimation 
reduces  and  increasing  the  phase  error  further  reduces  the  probability.  These 
results  indicated  that  phase  errors  need  to  be  minimized  as  much  as  possible 
during  the  hardware  design  as  well  as  during  DF  system  operation  to  improve  the 
accuracy  of  the  DF  system. 
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Figure  35.  Probability  of  correct  AOA  estimation  with  phase  errors. 

D.  EFFECTS  OF  FREQUENCY  CHANGE 

In  all  the  previous  simulations,  it  was  assumed  that  the  frequency  of  the 
signal  was  fixed  at  2.4  GHz.  In  order  to  investigate  into  how  sensitive  the 
baseline  design  was  to  changes  in  the  frequency,  the  physical  spacing  of  the 
antenna  elements  was  fixed  at  d- |  =30.93  cm,  d2  =10.44  cm  and  d3  =  3.61  cm. 

The  signal  frequency  was  varied  to  investigate  its  effect  on  the  performance  of 
the  baseline  design.  It  was  observed  in  Figure  36  that  tolerance  up  to  ±200 
MHz  bandwidth  is  possible  without  affecting  the  accuracy  of  the  transfer  function. 
As  the  signal  frequency  deviates  away  from  the  design  frequency,  the 
corresponding  transfer  function  also  deviates  accordingly.  This  can  be  seen  in 
Figure  37  where  the  slope  of  the  transfer  function  for  a  2.0  GHz  signal  no  longer 
resembles  the  slope  of  the  transfer  function  for  a  2.4  GHz  signal.  The  change  in 
slope  is  predictable  and  can  be  compensated  in  the  processing. 
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Figure  36.  Transfer  function  with  noise  using  2.4  GHz  and  2.2  GHz  signal. 


Figure  37.  Transfer  function  with  noise  using  2.4  GHz  and  2.0  GHz  signal. 


As  Figure  38  shows,  there  is  no  change  in  the  probability  of  correct  AOA 
estimation  even  with  a  2.0  GHz  signal.  To  design  for  a  wideband  frequency 
RSNS  DF  array  would  require  several  antenna  baselines  to  cover  the  desired 
bandwidth  with  the  necessary  accuracy. 
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Figure  38.  Probability  of  correct  AOA  estimation  baseline  design 

with  2.0  GHz  signal. 

E.  ALTERNATE  MODULI  SET 

The  baseline  design’s  low  value  moduli  set  was  selected  after  taking  into 
consideration  quantization  errors  in  a  low  SNR  scenario.  The  drawback  of  such 
a  low  moduli  set  was  a  small  dynamic  range,  which  resulted  in  lower  resolution. 
To  understand  the  effect  of  moduli  change,  an  alternate  higher  moduli  set  was 
chosen  for  comparison  while  keeping  the  same  h  ratio.  Table  15  shows  the 
alternate  moduli  set  and  some  of  the  key  parameters  obtained. 
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Moduli  set 

Dynamic 

Range 

Antenna  spacing 

FOV  (degree) 

[4  12  36] 

M Real  =113 

d- |  =  29.43  cm 
d2  =9.81  cm 
d3  =  3.27  cm 

±90 

[10  20  24] 

^Virtual  1  =  1 79 

d\  =  39.24  cm 
d’2=  19.62  cm 
c/3  =16.35  cm 

±28.4 

[20  40  48] 

^Virtual  2  = 

d\  =  39.24  cm 
d2  =19.62  cm 
do,  =16.35  cm 

±28.5 

Table  15.  Alternate  moduli  set’s  parameters. 


The  alternate  moduli  set  has  an  error  of  approximately  0.09  deg  RMS 
against  the  previous  baseline  design  error  of  0.5  deg  RMS.  The  increase  in 
accuracy  was  due  to  the  higher  dynamic  range  of  the  new  moduli  set.  Using  a 
similar  scenario,  the  transfer  function  and  the  probability  was  plotted  in  Figure  39 
and  Figure  40.  During  multiple  simulation  runs,  the  occurrence  of  gross  error 
was  observed  more  often  than  the  baseline  design.  This  observation  showed  up 
in  the  resulting  lower  probability  for  correct  AOA  estimation  in  Figure  40.  While 
the  higher  moduli  set  increases  resolution,  it  trades  off  on  the  overall  probability 
of  correctly  predicting  the  AOA. 
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Figure  39.  Transfer  function  of  alternate  moduli  set. 


Figure  40.  Probability  of  correct  AOA  estimation  for  alternate  moduli  set. 

The  effects  of  spacing  errors  and  phase  errors  were  also  plotted  for  the 
alternate  moduli  set  in  Figure  41  and  Figure  42.  The  tolerance  to  spacing  error 
increased  but  the  tolerance  to  phase  error  decreased.  Using  a  2.0  GHz  signal, 
the  alternate  moduli  set  exhibits  the  similar  effects  of  frequency  changes  as  the 
baseline  design.  This  can  be  seen  in  Figure  43. 
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Figure  41.  Transfer  function  for  alternate  moduli  with  spacing  errors. 


Figure  42.  Transfer  function  for  alternate  moduli  with  phase  errors. 
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Figure  43.  Transfer  function  for  alternate  moduli  using  2.4  GHz 

and  2.0  GHz  signal. 

F.  EFFECTS  OF  SNR  CHANGE 

If  the  demodulator  board  is  replaced  with  another  of  higher  sensitivity,  the 
baseline  design  would  be  able  to  handle  a  much  lower  SNR.  This  would  enable 
the  RSNS  virtual  DF  system  to  increase  its  maximum  range  of  detection.  To 
investigate  the  effects  of  such  low  SNR,  a  one  watt  emitter  is  simulated  at  30  km 
and  60  km.  The  probability  of  correct  AOA  estimation  was  plotted  for  both  the 
baseline  design  and  the  alternate  moduli  set  design.  While  both  designs  showed 
a  reduction  in  probability,  the  alternate  moduli  set’s  performance  at  longer  range, 
which  translates  to  low  SNR,  was  inferior  to  the  baseline  design.  This  can  be 
clearly  seen  by  comparing  Figure  44  and  Figure  45  for  a  signal  at  60  km.  The 
probability  of  correct  AOA  estimation  for  the  baseline  design  was  more  than  50% 
while  it  was  less  than  50%  for  the  alternate  moduli  set. 

In  order  to  select  the  correct  moduli  set,  the  most  probable  SNR  condition 
needs  to  be  assessed.  If  a  low  SNR  condition  is  expected,  the  baseline  design 
could  be  used.  If  a  high  SNR  condition  is  expected,  the  alternate  moduli  set 
should  be  used  instead. 
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Figure  44.  Probability  of  correct  AOA  estimation  for  baseline  design  with 

emitter  at  various  ranges. 


Figure  45.  Probability  of  correct  AOA  estimation  for  alternate  moduli  with 

emitter  at  various  ranges. 


G.  USING  TWO  SET  OF  MODULI 

Since  the  RSNS-based  interferometer  was  relatively  insensitive  to  small 
spacing  error,  an  investigation  into  using  two  complete  sets  of  different  moduli 
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based  on  the  same  ratio  was  carried  out.  Based  on  the  ratio,  h  =  y^,  the 

baseline  moduli  set  and  alternate  moduli  set  was  selected.  Using  the  baseline 
design’s  spacing  distance,  the  alternate  moduli  set  was  used  together  with  the 
baseline  moduli  set.  Although  this  effectively  introduced  spacing  error  into  the 
alternate  moduli  set,  the  limited  impact  on  the  accuracy  of  the  transfer  function 
may  be  worth  the  tradeoff  for  higher  probability  of  correct  estimation  in  low  SNR 
conditions.  The  transfer  function  for  the  baseline  and  alternate  moduli  set  is 
shown  in  Figure  46. 


Figure  46.  Transfer  function  for  baseline  and  alternate  modui  set. 
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VIII.  CONCLUSION 


This  thesis  focused  on  the  design  and  analysis  of  a  three-channel  RSNS 
virtual  spacing  DF  system.  By  analyzing  a  chosen  design  with  addition  of 
Gaussian  noise  as  well  as  carrying  out  Monte  Carlo  simulations,  the  sensitivities 
to  spacing  error,  phase  error,  frequency  error,  moduli  set  change  and  SNR 
change  were  determined.  The  baseline  design  showed  good  tolerance  to  most 
of  these  errors  and  should  be  considered  for  situations  where  low  SNR  is 
expected.  For  situations  where  high  SNR  is  expected,  the  alternate  moduli  set 
should  be  selected.  As  the  RSNS-based  DF  system  was  robust  to  the  small 
spacing  error,  the  best  solution  could  be  the  dual  use  of  both  low  and  high  moduli 
set.  This  would  give  a  high  probability  of  correct  estimation  at  all  time  while 
retaining  the  accuracy  of  correcting  mapping  of  the  estimated  AOA  to  the  actual 
AOA. 

The  application  of  RSNS  virtual  spacing  to  DF  has  been  demonstrated 
both  theoretically  as  well  as  experimentally  in  the  past.  It  showed  great  promise 
in  applying  this  concept  to  platforms  that  have  limited  size  but  yet  require  high 
accuracy  for  DF.  From  a  military  perspective,  an  unmanned  aerial  vehicle  used 
to  detect  a  low  probability  of  intercept  emitters  is  one  of  the  likely  scenarios 
where  RSNS-based  DF  concept  could  be  applied. 

Some  of  the  possible  future  research  areas  to  continue  this  RSNS-based 
DF  concept  are: 

1.  To  build  the  baseline  design  using  available  hardware  and  verify 
the  performance  of  the  design  with  its  theoretical  performance.  With 
the  experimental  data,  other  errors  could  be  identified  and 
corrected.  Initial  testing  could  be  done  in  the  laboratory  before 
expanding  it  to  field  testing. 

2.  To  determine  the  best  set  of  two  or  more  moduli  set  in  order  to 
optimize  the  performance  of  the  RSNS-based  DF  system. 
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3.  Further  improvement  of  the  simulation  by  including  the  effects  of 
mutual  coupling  interference. 

4.  The  use  of  non-PRP  moduli  set  resulted  in  a  smaller  dynamic  range 
than  what  would  have  been  possible  with  a  PRP  moduli  set. 
Studies  into  the  use  of  PRP  moduli  set  for  virtual  spacing  should 
continue  in  order  to  fully  exploit  the  advantage  RSNS-based  DF 
system. 
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